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Events containing hadronically decaying heavy particles with large momentum, leading to so- 
called merged jets, are expected to play a significant role in both searches for new physics and 
measurements of Standard Model processes at the CERN Large Hadron Collider. In this article, 
a comparative study of the modeling of such merged jet topologies by different Monte-Carlo event 
generators is presented. The observed differences emphasize the need to refine such modeling based 
on the observation of Standard Model processes prior to claims of discovery. 



I. BACKGROUND 

The increased energy and luminosity of modern hadron 
colliders, first at Run II of the Fermilab Tevatron 
Collider and now at the CERN Large Hadron Collider 
(LHC), requires increasingly sophisticated Monte Carlo 
(MC) generators to model both background and signal 
processes in order to extract new physics results from the 
data. Indeed, one major challenge facing these analyses 
is that of large final state jet multiplicities, with events 
containing 6 or more high transverse momentum (jpt) 
jets expected to be common at the LHC. Additionally, 
jet structure is likely to play an increasingly important 
role in physics analyses of hadronic decays of boosted 
massive particles [THS]. The latter aspect in particular 
probes a difficult region of phase space, testing the ability 
of MC generators to accurately simulate jet production. 

Monte Carlo event generation is commonly performed 
in two main stages. A "matrix element step" evaluates 
a hard process with a fixed number of incoming and 
outgoing particles (e.g. 2^2), while additional jets 
may be produced during the parton shower phase of 
the generator. While it would be infeasible to evaluate 
the entire event using the hard process matrix element 
methods, it is generally difficult to accurately generate 
multiple additional hard jets in the parton shower phase. 

In order to improve the modeling of events with 
multiple hard jets, a number of next generation MC 
generators have been developed to generate hard pro- 
cesses with an increasing number of outgoing particles 
(up to 2 ^ 9) during the matrix element step. However, 
the same jet may now be produced in either the matrix 
element or parton shower, creating significant ambiguity 
in the division of phase space for jet production and 
the potential for "double counting." Thus, a proper jet 
matching algorithm which can dictate which jets should 
be produced in the hard process and which in the parton 
shower is crucial to ensure proper phase space coverage, 
a task which is made difficult by the intrinsic difficulty in 
linking hard process quarks or gluons to final- state jets. 

This article compares the performance of several MC 
generators. Pythia[7, 8 and Herwig++[9 are tradi- 



tional MC generators which generate only 2^2 hard 
processes and require any additional jet production to 
come from the parton shower. Alpgen[TO], a 2 ^ n 
in < 9) generator implements MLM matching, which 
allows the event evolution to proceed without restriction 
but afterwards vetoes events whose hard jets do not 
match the parton-level quarks and gluons produced in 
the hard process. CKKW[11 matching, which suppresses 
the production of soft jets in the hard process and hard 
jets during the parton shower phase according to the 
scale of each individual branch splitting as the event is 
being generated, is implemented in Sherpa[12 , also a 
2 ^ n (n < 9) generator. 

Generators such as Alpgen and Sherpa offer addi- 
tional advantages when considering the decays of heavy 
particles. While 2^2 generators can generate heavy 
particles in the matrix element, heavy particle decays are 
usually handled at a later phase in the event evolution, 
before the parton shower. However, such a factorization 
is only approximate, and modeling may be improved by 
incorporating heavy particle decays into the hard process, 
as is possible in Alpgen and Sherpa. Decay products 
from heavy particles are not subject to double counting 
and are generally excluded from the jet matching proce- 
dure. 

Hard jets produced in the parton shower or from 
radiation from the hard process, rather than from the 
decay of a heavy particle, are hereafter referred to as 
QCD jets. In either Alpgen or Sherpa, the user must 
specify a maximum number of such QCD jets which may 
be produced in the hard process. For events at this upper 
limit, additional jet production in the parton shower 
is permitted provided that it is softer than the softest 
jet from the hard process (according to an algorithm- 
dependent metric.) Thus, both of these approaches are 
designed to produce full phase space coverage without 
double counting. 

Previous studies of these and other generators have 
compared the jet pt and rapidity distributions in pp and 
pp collisions [13 as well as distance [14| in e^e~ 
jets events. This study utilizes pp ^ ti ^ jets events 
to examine the structure created through the decay of 
boosted heavy particles, while high multiplicity pp 
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jets (QCD) events provide a more general source of jet 
mergers. Both data sets are sensitive to characteristics 
of the individual matrix element (ME) generators and 
the parton showers, while the second set in particular 
may probe differences between CKKW and MLM jet 
matching. 



II. JET MERGING AND MATCHING 
ALGORITHMS 

In events with a large jet multiplicity, it will sometimes 
occur that a pair of 'low mass jets' has little separation. 
This may result either from the decay kinematics of a 
boosted heavy particle or from a large multiplicity of 
uncorrelated jets. In such cases, jet reconstruction may 
merge this pair of 'low mass jets' into a single jet which 
acquires a large mass. This implicitly leads to the appear- 
ance of a peak or shoulder in the jet mass distribution 
whose properties depend on the jet algorithm and its 
parameters, characteristics of the parton shower, and the 
kinematics of any heavy particle decays. Thus, the jet 
mass distribution provides an experimentally accessible, 
sensitive probe of the effects of jet merging and the 
performance of MC generation tools for modeling jet 
structure. 

In MLM matching, particle-level jets reconstructed 
with a simple cone algorithm are matched to hard quarks 
or gluons from the matrix element phase. In the case 
that two hard QCD jets are merged by this algorithm, 
a mismatch between the observed and expected jets will 
trigger a veto unless exactly one of the jets originates 
from the ME hard process and the other the parton 
shower. On the other hand, hard jets radiated during 
the parton shower may escape the veto provided they are 
sufficiently close to a hard process QCD jet so as to be- 
come merged (according to the internal cone algorithm.) 
However, mapping of hard quarks or gluons to particle- 
level jets has limited accuracy due to QCD confinement 
and higher order effects, potentially resulting in improper 
event vetoes. 

Conversely, in CKKW two jets from the hard process 
may be allowed to merge without causing a veto. How- 
ever, every jet radiated during parton showering must be 
soft in relation to its point of origin even if it merges with 
another jet. Thus, merged hard jets may only originate 
from the hard process. 

In either case, merged jets may be considered to exist 
along a boundary in phase space between an N and 
N-hl jet hard process. Therefore, jet merging occurs in 
a region of phase space which is particularly sensitive 
to the phase space coverage produced by jet matching, 
and the size and position of the shoulder in the jet 
mass distribution is an excellent probe of the differences 
between jet matching algorithms (and reality.) 
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FIG. 1: W boson pt distribution in ti events for 
several generators. 

III. tt EVENT GENERATION AND JET 
RECONSTRUCTION 

For these studies pp tt^ jets events were generated 
with several generators at a 14 TeV center-of-mass energy 
and with the requirement that each top quark decay 
hadronically into 3 jets via a W boson. Jet multiplicities 
of 7-8 were thus achieved, providing a jetty environment 
in which mergers become likely. No generator-level cuts 
were applied and multiple interactions were turned off in 
every generator. 

The following generators were used: 

• Pythia 6.325, using the CTEQ5L parton distribu- 
tion functions (PDFs.) 

• Pythia 8.150, using the CTEQ5L PDFs. 

• HERWIG++ 2.4.2, using the 2008 MRST leading 
order PDFs. 

• Alpgen 2.13 for ME generation using the CTEQ5L 
PDFs and Pythia 6.325 for parton showering and 
hadronization. Only tt ^ 1 (exclusive) or 2 (in- 
clusive) light parton processes were considered, as 
the production cross section for + light parton 
events is vanishingly small. MLM matching was 
applied to ensure proper phase space coverage, and 
the two samples were manually merged according 
to the cross-section estimates provided by Alpgen. 
By default, Alpgen generates only on-shell top 
quarks and W bosons. 

• Sherpa 1.2.2 using the COMIX ME generator, 
which is recommended for large particle multiplici- 
ties, and the CTEQ6L PDFs. The CKKW merging 
scale was set to 30 GeV. tt + 0, 1, and 2 light 
parton hard processes were considered and com- 
bined internally by the Sherpa generator. Parton 
showering and hadronization were also performed 
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FIG. 



2: pt distribution of all jets with \y\ < 4.5, pr > 70 GeV for several jet clustering algorithms. Left to right: 
DO Run II Cone R=1.0, anti-A:^ R=0.8, Cambridge/ Aachen R=1.0 
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3: Pt distribution of merged W^-jets with \y\ < 4.5, pr > 70 GeV for several jet algorithms. Left to right: DO 
Run II Cone R=1.0, anti-/c^ R=0.8, Cambridge /Aachen R=1.0 



by Sherpa. Top quarks and W bosons were forced 
to be on shell. 

Jets were reconstructed at particle-level (i.e. after 
parton shower and hadronization) using the DO Run 
II Cone [15^, anti-/cx[16^, and Cambridge/Aachen[17 jet 
algorithms with multiple radii. This was performed using 
Spartyjet 3.4.1[TF and Fastjet 2.4.1[T9^. Addition- 
ally, parton-level information about the hard process 
W bosons and their immediate decay products were 
extracted from the event history. 



IV. tt ANALYSIS 

Figure [l] shows the particle-level W boson pt distribu- 
tions for each generator, with generally good agreement 
but some divergence at higher pt- However, Fig. |2] 
shows differences in the jet pr spectrum for all jets with 
\y\ < 4.5, Pt > 70 GeV. 



A. Merged W-Jets 

Jets containing both products from a boosted hadronic 
W boson decay (hereafter referred to as 'merged W-jets') 



are considered separately in this study due to their value 
as an experimental calibration signal for future studies 
on jet structure. Additionally, the strong kinematic 
constraints on such jets generate distinctive features 
in their mass distribution which must be considered 
separately from other jet mergers. Such jets are tagged 
via the requirement that each of the parton-level decay 
products from the W boson be sufficiently close to the jet 
axis (v/At/2 + A(/)2 < R.) The remaining bottom quark 
from the top quark decay is allowed but not required to 
be contained in the jet as well. 

Figure |3] shows a non-trivial difference in the the 
Pt distribution for merged W-jets. The jet mass dis- 
tributions shown in Fig. [4] display a pair of peaks at 
approximately 80 and 170 GeV, corresponding to the 
masses of the W boson and top quark respectively. 

There are several differences between the MC genera- 
tors which may be seen in these distributions, especially 
in the treatment of the W boson mass spectra. Sherpa 
and Alpgen force heavy particle decays (both W bosons 
and top quarks) to be on the mass shell (necessary to keep 
the computational cost manageable given the increased 
responsibilities of the hard process.) This results in a 
sharper W boson mass peak in Sherpa and Alpgen 
than is evident in the 2^2 generators. 

On the other hand, both HERWIG++ and Pythia6 
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FIG. 4: Jet mass distribution of merged W-jets with \y\ < 4.5^ pr > 70 GeV for several jet algorithms. Left to 
right: DO Run II Cone R=1.0, anti-A:^ R=0.8, Cambridge/ Aachen R=1.0 



generate the W boson according to a relativistic Breit- 
Wigner mass distribution with a fixed width, while 
PythiaS implements a running width. This results in 
a sharp suppression of the low mass tail in PythiaS 
compared to the other 2^2 generators, and a larger tail 
in the high mass region. Although only a fraction of the 
total, W bosons at the low end of the mass spectrum are 
more likely to produce merged jets than W bosons which 
are near or above the mass shell, especially at moderate 
Pt. This results in a significant excess of 'merged W jets' 
in the low mass low pt turn on regions in HerwigH — h 
and PythiaG when compared to the other generators. 
It is not fully understood whether it is better to use a 
fixed or running mass width in these decays, and these 
differences should be viewed as a systematic uncertainty 
in the modeling of top quark decays and other heavy 
particles [20 . 

However, differences in particle width modeling do not 



appear to have as strong of an effect on the shape of 
the second mass peak, as Sherpa and Alpgen do not 
show narrower distributions than the other generators 
despite generating top quarks on the mass shell. Rather, 
in this case Pythia6, Herwig++ and Alpgen, ah of 
which use virtuality ordered parton showers, show close 
agreement in the shape of this peak when compared 
to PythiaS and Sherpa which shower in pt and kj_ 
respectively. 

It is important to note that the generators may show 
differences in the documentation of intermediary parti- 
cles in the event history which may infiuence the tagging 
of merged W-jets. Since such intermediary particles 
are inherently unphysical, these differences should not 
be refiected in the physics predictions produced by the 
generators. However, it is possible that they would 
infiuence the tagging criteria in this analysis. Figure [5] 
shows the effect of tightening the selection criteria on 
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FIG. 5: Jet pt and mass distribution of jets tightly pointed at by a boson with y < 4.5, pr > 70 GeV. Both W 
boson daughter quarks are required to be within R/2 from jet axis. Jets are reconstructed with Cambridge /Aachen 

R=1.0 
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FIG. 6: Jet mass distribution of merged VF-jets with \y\ < 4.5, pr > 70 GeV in Pythia6, PythiaS (a) and 
Alpgen (b) before and after reweighting jets to match PythiaG's pt distribution, using Cambridge /Aachen R=1.0 

jets. 



merged W-jets by requiring both quarks to be within 
half the jet radius. The relative insensitivity of these 
distributions to the tagging criteria suggests that the 
observed disagreements reflect differences in the final 
states produced by the generators, rather than internal 
differences in the event record. This tightening does 
however suppress the low mass, low pr end of the 
distributions, which in turn makes the influence of the 
W boson width modeling more obvious. 

Differences in the jet pt spectrum between generators 
account for only a small part of the differences in the 
mass spectrum. Reweighting the merged 1^-jets such 
that the jet pr distributions match does lead to some 



improve the agreement between generators but does not 
fundamentally alter the behavior, as shown in Fig. |6] 

Increasing the jet radius suppresses the first peak in 
the mass spectrum and fills in the gap between the two 
peaks, suggesting that it becomes more difficult to create 
a 'merged W^-jet' which does not include at least some 
additional radiation, perhaps from the third jet. This 
is shown in Fig. [7] for Pythia 8, Alpgen and Sherpa. 
Differences in behavior between the generators as the jet 
radius is varied would suggest that different generators 
may lead to different jet radius optimizations for the 
tagging of heavy particle decays. 
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FIG. 7: Jet mass distribution of merged W-jets with \y\ < 4.5, pr > 70 GeV in PythiaS, Alpgen and Sherpa for 

several jet radii. 
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FIG. 8: Mass distribution of jets not pointed at by a boson with \y\ < 4.5, pr > 70 GeV for several jet 
clustering algorithms. Left to right: DO Run II Cone R=1.0, anti-/cx R=1.0, Cambridge /Aachen R=1.0 
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FIG. 9: Mass distribution of jets not pointed at by a 
W boson with \y\ < A.b.pr > 70 GeV for the anti-/c^ 
R=0.8 jet algorithm after jet pt reweighting. 




FIG. 10: Mass distribution of jets not pointed at by a 
W boson with \y\ < 4.5, > 70 GeV for the anti-/c^ 
R=0.8 jet algorithm. Merged b-jets are those which are 
pointed at by a bottom quark and exactly one W boson 

decay product from the same top quark. Non-merged 
jets include all jets which are neither merged W-jets nor 
merged b-jets (but may still contain jet mergers). 
Shown for HERWIG++. 
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FIG. 11: pt and mass distribution of the leading jet {\y\ < 1.5, > 300 GeV) reconstructed with the anti-/cx, 
R=1.0 jet algorithm for QCD multijet events with multiple 'narrow' jets in a tight range \y\ < 1.5. Narrow jet 
preselection is based on either anti-/cx (a,c) or DO Run II Cone (b,d) with R=0.2. 



B. Bottom Jet Mergers 

Merged jets which contain the majority of the radiation 
produced in the decay of a single heavy particle (such as 
a top quark or W boson) are subject to tight kinematic 
constraints. However, it is also possible that the bottom 
quark resulting from a top quark decay will merge with 
exactly one of the W boson decay products from the 
same top quark (hereafter referred to as 'merged 6-jets.') 
Such mergers are less constrained by decay kinematics, 
limiting their usefulness as a calibration tool. However, 
such mergers may prove to be a non-trivial background 
in other physics analyses. 

Figure [S] shows the mass distributions for all jets 
excluding merged W^-jets and requiring pt > 70 GeV, 
\y\ < 4.5 for each generator and several jet algorithms. 



The small bump in the distribution produced by the 
DO Run II Cone algorithm at around 170 GeV suggests 
some inefficiency in the tagging of such mergers, possibly 
because this algorithm absorbs nearby radiation easily, 
producing a larger effective jet area for a given radius. 
Additionally, the cone algorithm appears to be less 
sensitive to differences between the generators in this and 
several previous distributions. 

There is a small shoulder in these mass distributions 
appearing between 100 and 120 GeV, especially for the 
anti-/c^ algorithm. This shoulder appears to share a 
similar position for each generator, and is most apparent 
in those generators with the softest jet mass spectra. The 
disagreement between generators in these distributions 
is at least partially explained by differences in the pt 
spectrum. This is demonstrated in Fig. [9j which shows 
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FIG. 12: Mass distribution of the leading jet 
(\y\ < 1.5, > 300 GeV) reconstructed with the 
anti-Zc^, R=1.0 jet algorithm for QCD events with 
multiple 'narrow' (DO Run II Cone R=.2) jets. Shown 
for the Sherpa MC generator, before and after 
reweighting by leading jet pt- 



the mass distributions for non-' merged W jets' after these 
jets are reweighted such that the jet pr distribution 
matches PythiaG. However, after such a reweighting is 
applied, the mass shoulder appears to be more significant 
in those generators with the softest original jet mass 
spectrum, especially HerwigH — h- 
Figure 



10 shows the effect of separating out 'merged 



b-jets' from the entire non-'merged W jet' distribution. 
This shows that the shoulder in the mass distribution is 
primarily produced by mergers between a bottom quark 
and the W boson decay product from the same top quark. 



V. QCD EVENT GENERATION AND 
RECONSTRUCTION 

The preceding sections discussed the production of 
merged jets which are produced in the decay of boosted 
heavy particles. The following two sections will explore 
more general jet mergers which may be produced in any 
event with a large jet multiplicity. 

For this purpose, QCD multijet events were generated 
using the same set of generators as described in section 
III (Pythia 8.135 was used instead of 8.150.) Only weak 
generator-level cuts were applied in order to minimize the 
introduction of bias, as cuts available to 2 ^ 2 generators 
are generally not equivalent to those available for 2 ^ n 
processes. In each 2^2 generator, pt > 100 GeV was 
required in the hard process. The following generator 
settings were used in Alpgen and Sherpa: 

• Alpgen: 2-5 jet hard processes were generated 
with the requirement that all partons have Pt > 
40 GeV, \y\ < 2.5. After parton showering, MLM 



matching was applied at a scale of 48 GeV and 
requiring jet \y\ < 2.0. 

• Sherpa: 2-5 jet hard processes were generated 
using the COMIX ME generator, with the CKKW 
merging scale at 30 GeV. The two primary parton- 
level jets in the event are required to have pt > 
60 GeV. The leading particle- level anti-Zcx = 1.3 
jet in the event with \y\ < 2.0 is required to have 
Pt > 200 GeV. 

As in section III, multiple interactions were turned 
off in every generator and jets were reconstructed using 
several jet algorithms. 



VI. QCD ANALYSIS 

This analysis of QCD events proceeds by selecting 
events with a large multiplicity of 'narrow' (small radius) 
jets before applying a 'fat' (larger radius) jet algorithm 
and observing the resulting mass distribution. The 
performance of the MC generators is then evaluated by 
observing the degree of agreement between them. In this 
section, two variations of this analyses will be presented. 

In the first variation, events are preselected if they 
contain at least four 'narrow' jets (ie. jets which are 
reconstructed with either the DO Run II Cone or anti-/cx 
algorithm with R=.2) with \y\ < 1.5 and pt > 60 GeV, 
and at least one with pt > 170 GeV. Events are then 
required to have a leading (anti-/cx with R=1.0) jet with 
Pt > 300 GeV and \y\ < 1.5. 



Figure 11 shows the leading jet pt and mass distri- 
butions of events surviving this selection. The mass 
distribution produced by the DO Run II Cone narrow 
jet preselection shows a clear shoulder produced by the 
presence of merged narrow jets, while the distribution 
produced by anti-/cx preselection suggests a similar, 
although less pronounced effect. 

Pythia8 appears to show some disagreement with the 
other generators, with a slightly slower turn on under 
either choice of preselection 'narrow' jet definition. This 
has the effect of reducing the apparent dip between 
the two mass peaks when the DO Run II Cone narrow 
jet algorithm is used, and creating a more apparent 
shoulder when the anti-/cx narrow jet preselection is used. 
Sherpa appears to have a harder spectrum than the 
other generators, and Fig. [12] shows that reweighting 
Sherpa's Pt spectrum to match HerwigH — h improves 
the agreement in the mass distribution. 

A second analysis technique leads to a more direct 
selection of merged jets. This is achieved by selecting 
a leading fat (anti-/c^ with R=1.0) jet with \y\ < 1.5 
and Pt > 300 GeV and requiring at least 2 narrow jets 
with Pt > 60 GeV be contained within its jet radius. 



The result of this approach can be seen in Fig. 13, which 
shows the leading jet pt and mass distributions for both 
narrow jet definitions. 
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FIG. 13: and mass distribution of the leading jet {pt > 300 GeV, \y\ < 1.5) in QCD multijet events, 
reconstructed with the anti-/cx, R=1.0 jet algorithm and requiring that multiple 'narrow' jets be contained within 
the jet radius. Narrow jet preselection is based on either anti-/cx (a,c) or DO Run II Cone (b,d), R=0.2. 



Such 'merged QCD' jets form a nontrivial fraction of 
all leading QCD jets. Table |T] displays the efficiency 
for events in each generator which contain at least one 
fat jet with pt > 300 GeV and \y\ < 1.5 to survive 
these selection cuts. Furthermore, although the shapes 
of the distributions shown in Fig. [13] are similar in each 
generator, the efficiencies vary significantly, highlighting 
the need to better understand the production of such 
QCD jets as a potential background in other analyses 
utilizing jet structure. 



VII. CONCLUSIONS 

This study has examined the modeling of merged 
jets both from boosted particle decays and from QCD 



multijet backgrounds. Such jet mergers provide a window 
into MC jet modeling and are a sensitive probe of jet 
matching, parton showering, and other characteristics 
of MC generators. While decent overall agreement 
exists between the generators, significant differences, 
both explained and unexplained, were observed in key 
observables such as jet mass, with substantial variations 
in event rates observed in several regions of phase 
space. Such differences may provide a handle towards 
a better understanding of systematic uncertainties which 
may affect precision measurements and searches for new 
physics at the LHC, especially those involving the decay 
of boosted heavy particles, and may prove to be a 
useful starting point for future studies to improve the 
performance of modern MC generators by comparing 
these predictions to LHC data as it becomes available. 
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Selection 1 


Selection 2 


narrow jet : 


Anti-A:± 


DO Run II Cone 


Anti-A:± 


DO Run II Cone 


Pythia6 

HERWIG + + 

Alpgen 
PythiaS 
Sherpa 


.137 ± .002 
.126 =b .002 
.133 ± .003 
.163 ± .002 
.181 ± .002 


.070 ± .001 
.066 ± .002 
.071 ± .001 
.128 ± .002 
.109 ±.002 


.048 ± .001 
.051 =b .001 
.045 ± .002 
.078 ± .001 
.057 zb .001 


.019 ± .001 
.021 dz .001 
.015 zb .001 
.038 dz .001 
.028 ± .001 



TABLE I: Efficiency for events with at least one jet with pr > 300 GeV and \y\ < 1.5 to survive event preselection 

criteria for each generator. 
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